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ABSTRACT 

In this lecture I give an overview of shock acceleration, interactions of high energy 
cosmic rays with, and propagation through, the background radiation, and the 
resulting electron-photon cascade. I argue that while the origin of the highest 
energy cosmic rays is still uncertain, it is not necessary to invoke exotic models 
such as emission by topological defects to explain the existing data. It seems likely 
that shock acceleration at Fanaroff-Riley Class II radio galaxies can account for 
the existing data. However, new cosmic ray data, as well as better estimates of the 
extragalactic radiation fields and magnetic fields will be necessary before we will 
be certain of the origin of the highest energy particles occurring in nature. 



1 Introduction 

Cosmic rays with energies up to 100 TeV are thought to arise predominantly 
through shock acceleration by supernova remnants (SNR) in our Galaxy 
A fraction of the cosmic rays accelerated should interact within the supernova 
remnant and produce gamma-rays ^, and recent observations above 100 
MeV by the EGRET instrument on the Compton Gamma Ray Observatory 
have found gamma ray signals associated with at least two supernova remnants 
- IC 443 and 7 Cygni (however, it is possible that the gamma ray emission 
from IC 443 is associated with a pulsar within the remnant rather than the 
remnant itself |^). Further evidence for acceleration in SNR comes from the 
recent ASCA observation of non-thermal X-ray emission from SN 1006 
Reynolds and Mastichiadis Q interpret the latter as synchrotron emission 
by electrons accelerated in the remnant up to energies as high as 100 TeV, 
although Donau and Biermann |9| suggest it may be bremsstrahlung from 
much lower energy electrons. 

Acceleration to somewhat higher energies than 100 TeV may be possible 
[01, but probably not high enough to explain the smooth extension of the 
spectrum to 1 EeV. Several explanations for the origin of the cosmic rays 
in this energy range have been suggested: reacceleration of the supernova 
component while still inside the remnant [Tl|] ; by several supernovae exploding 
into a region evacuated by a pre-supernova star ||l^ ; or acceleration in shocks 
inside the strong winds from hot stars or groups of hot stars ||l^ . At 5 EeV the 
spectral slope changes, and there is evidence for a lightening in composition [|l^ 
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and it is likely this marks a change from galactic cosmic rays to extragalactic 
cosmic rays being dominant. 

The cosmic ray air shower events with the highest energies so far detected 
have energies of 2 x lO" GeV and 3 x lO" GeV Q. The question of 
the origin of these cosmic rays having energy significantly above 10"'^^ GeV is 
complicated by propagation of such energetic particles through the universe. 
Nucleons interact with the cosmic background radiation fields, losing energy 
by pion photoproduction, and may emerge as either protons or neutrons with 
reduced energy. The threshold for pion photoproduction on the microwave 
background is ~ 2 x 10^° GeV, and at 3 x 10^^ GeV the energy-loss distance is 
about 20 Mpc. Propagation of cosmic rays over substantially larger distances 
gives rise to a cut-off in the spectrum at ^ 10^^ GeV as was first shown 
by Greisen and Zatsepin and Kuz'min ]T^, the "GZK cut-off", and a 
corresponding pile- up at slightly lower energy p9|, [2^ . These processes occur 
not only during propagation, but also during acceleration and may actually 
limit the maximum energies particles can achieve. 

In this lecture I give an overview of shock acceleration, describe interac- 
tions of high energy protons and nuclei with radiation, discuss maximum ener- 
gies obtainable during acceleration, outline propagation of cosmic rays through 
the background radiation and the consequent electron-photon cascading, and 
finally discuss conventional and exotic models of the highest energy cosmic 
rays. 



2 Cosmic Ray Acceleration 

For stochastic particle acceleration by electric fields induced by motion of mag- 
netic fields B, the rate of energy gain by relativistic particles of charge Ze can 
be written (in SI units) 

§ - ^Zec^B (1) 

acc 

where ^ < 1 and depends on the acceleration mechanism. I shall give a sim- 
ple heuristic treatment of Fermi acceleration based on that given in Gaisser's 
excellent book [^. I shall start with 2nd order Fermi acceleration (Fermi's 
original theory) and describe how this can be modified in the context of super- 
nova shocks, or other strong astrophysical shocks, into the more efficient 1st 
order Fermi mechanism at supernova (SN) or other shocks. More detailed and 
rigorous treatments are given in several review articles [22|-|2^] . See the review 
by Jones and Ellison on the plasma physics of shock acceleration which 
also includes a brief historical review and refers to early work. 
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2.1 Fermi's Original Theory 



Gas clouds in the interstellar medium have random velocities of ^ 15 km/s 
superimposed on their regular motion around the galaxy. Cosmic rays gain 
energy on average when scattering off these magnetized clouds. A cosmic ray 
enters a cloud and scatters off irregularities in the magnetic field which is tied 
to the cloud because it is partly ionized. 

In the frame of the cloud: (a) there is no change in energy because the 
scattering is collisionless, and so there is elastic scattering between the cosmic 
ray and the cloud as a whole which is much more massive than the cosmic ray; 
(b) the cosmic ray's direction is randomized by the scattering and it emerges 
from the cloud in a random direction. 




Figure 1: Interaction of cosmic ray of energy Ei with "cloud" moving with speed V 



Consider a cosmic ray entering a cloud with energy Ei and momentum pi 
travelling in a direction making angle 0i with the cloud's direction. After scat- 
tering inside the cloud, it emerges with energy E2 and momentum p2 at angle 
62 to the cloud's direction (Fig. |l|). The energy change is obtained by applying 
the Lorentz transformations between the laboratory frame (unprimed) and the 
cloud frame (primed). Transforming to the cloud frame: 

=7£;i (1-/3 cos 6*1) (2) 

where (3 = V/c and 7 = 1/ \/l — /3^. 
Transforming to the laboratory frame: 

E2=^E'^{1^ f3cose'^). (3) 

Since E'2 — E[ we obtain the fractional change in energy {E2 ~ Ei)/Ei, 

AE _ 1- I3COS01+ P cos 9'^ - cos 6*1 cos 6>^ 
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We need to obtain average values of cos 9i and cos Inside the cloud, the 
cosmic ray scatters off magnetic irregularities many times so that its direction 
is randomized, 

(cos6l^)=0. (5) 

The average value of cos^i depends on the rate at which cosmic rays collide 
with clouds at different angles. The rate of collision is proportional to the 
relative velocity between the cloud and the particle so that the probability per 
unit solid angle of having a collision at angle 0i is proportional to {v — V cos 9i). 
Hence, for ultrarelativistic particles {v = c) 

dP 

— cx (1-/3 cos (6) 



and we obtain 



givmg 



'^o.)-Jco.o,§-dn,/J§-dn, = -l, (7) 



1 ~ -(3' (8) 



E 1-^2 - 3^ 
since (3-^1. 

We see that {AE)/E cx (3"^ is positive (energy gain), but is 2nd order in 
f3 and because /3 <C 1 the average energy gain is very small. This is because 
there are almost as many overtaking collisions (energy lost) as there are head- 
on collisions (energy gain). 



2.2 1st Order Fermi Acceleration at SN or Other Shocks 

Fermi's original theory was modified in the 1970's p6|-p9t to describe more 
efficient acceleration (1st order in /?) taking place at supernova shocks but is 
generally applicable to strong shocks in other astrophysical contexts. 

During a supernova explosion several solar masses of material are ejected 
at a speed of ~ 10^ km/s which is much faster than the speed of sound in 
the interstellar medium (ISM) which is ^ 10 km/s. A strong shock wave 
propagates radially out as the ISM and its associated magnetic field piles up 
in front of the supernova ejecta. The velocity of the shock, Vs, depends on the 
velocity of the ejecta, Vp, and on the ratio of specific heats, 7. The SN will 
have ionized the surrounding gas which will therefore be monatomic (7 = 5/3), 
and theory of shock hydrodynamics shows that for a strong shock 

Vs/Vp ~ 4/3. (9) 
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In order to work out the energy gain per shock crossing, we can visuahze 
magnetic irregularities on either side of the shock as clouds of magnetized 
plasma of Fermi's original theory (Fig. |^). By considering the rate at which 
cosmic rays cross the shock from downstream to upstream, and upstream to 
downstream, one finds (cos 6*1) ~ —2/3 and (cos 6*2) = 2/3, giving 



(Ai?> 



4 Vs 
E 3^ c 



(10) 



Note this is 1st order in (3 and is therefore more efficient than Fermi's original 
theory. This is because of the converging flow - whichever side of the shock 
you are on, if you are moving with the plasma, the plasma on the other side 
of the shock is approaching you at speed Vp. 
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Figure 2: Interaction of cosmic ray of energy Ei with a shock moving with speed Vs. 

To obtain the energy spectrum we need to find the probability of a cos- 
mic ray encountering the shock once, twice, three times, etc. If we look at 
the diffusion of a cosmic ray as seen in the rest frame of the shock (Fig. 
there is clearly a net flow of the energetic particle population in the down- 
stream direction. The net flow rate gives the rate at which cosmic rays are lost 
downstream 

-Rioss = "cr1^s/4 ni-2s-i (11) 

since cosmic rays with number density ncR a-t the shock are advected down- 
stream with speed Vs/4: (from right to left in Fig. ||). 

Upstream of the shock, cosmic rays travelling at speed v at angle 6 to 
the shock normal (as seen in the laboratory frame) approach the shock with 
speed (Vs -|- t^cos^) as seen in the shock frame. Clearly, to cross the shock, 
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«2= Vs /4 u^= Vs 

downstream upstream 



Figure 3: Diffusion of cosmic rays from upstream to downstream seen in the shock frame 
speed Vs. 



COS 9 > —Vs/v. Then, assuming cosmic rays upstream are isotropic, the rate 
at which they cross from upstream to downstream is 

-Rcross = '^CRi- / {Vs + vcose)2Tvd{cos0) =ncRv/4: m ^ (12) 

47r j-vs/v 

The probabihty of crossing the shock once and then escaping from the 
shock (being lost downstream) is the ratio of these two rates: 

Prob. (escape) = -RiossZ-Rcross = Vs/v (13) 

where we have neglected relativistic transformations of the rates because Vs <^ 
c. The probability of returning to the shock after crossing from upstream to 
downstream is 

Prob. (return) = 1 — Prob. (escape) (14) 

and so the probability of returning to the shock k times and also of crossing 
the shock at least k times is 

Prob.(cross > fc) = [1 - Prob. (escape)]''. (15) 

Hence, the energy after k shock crossings is 

E = Eoi^l + —j (16) 

where Eq is the initial energy. 

To derive the spectrum, we note that the integral energy spectrum (number 
of particles with energy greater than E) on acceleration must be 

Q{> E)(x[l- Prob. (escape)] (17) 
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where 

"'^ \n{l + ^E/Ey 

Hence, 

lnQ(> E) = A+^^^^L-^ ln[l - Prob. (escape)], (19) 
where A is a constant, and so 

hiQ{> E) = B -V\tlE (20) 

where _B is a constant and 

^ ln[l- Prob. (escape)] _ 

in(i + A£:/i;) ■ ^ 

Hence we arrive at the spectrum of cosmic rays on acceleration 

Q(> E) cx E-'^ (integral form) (22) 

Q{E)(xE~^ (differential form). (23) 

The observed cosmic ray spectrum is steepened by energy- dependent escape 
of cosmic rays from the Galaxy. 

2.3 Shock Acceleration Rate 

The rate of gain of energy is given by 



dE 



^ (24) 

''cycle 



where tcycic is the time for one complete cycle, i.e. from crossing the shock 
from upstream to downstream, diffusing back towards the shock and crossing 
from downstream to upstream, and finally returning to the shock. We shall 
discuss this process in the shock frame (see Fig. ||) and consider first particles 
crossing the shock from upstream to downstream and diffusing back to the 
shock, i.e. we shall work out the average time spent downstream. Since we 
are considering non-relativistic shocks, the time scales are approximately the 
same in the upstream and downstream plasma frames, and so in this section I 
shall drop the use of subscripts indicating the frame of reference. 

Diffusion takes place in the presence of advection at speed U2 in the down- 
stream direction. The typical distance a particle diffuses in time t is ^/k2t 
where k2 is the diffusion coefficient in the downstream region. The distance 
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advected in this time is simply U2t. If ^/k2i U2t the particle has a very high 
probability of returning to the shock, and if \/k2t <C the particle has a very 
high probability of never returning to the shock (i.e. it has effectively escaped 
downstream). So, we set \/k2t = U2t to define a distance k2/u2 downstream 
of the shock which is effectively a boundary between the region closer to the 
shock where the particles will usually return to the shock and the region far- 
ther from the shock in which the particles will usually be advected downstream 
never to return. There are ncRk2/u2 particles per imit area of shock between 
the shock and this boundary. Dividing this by i?cross we obtain the average 
time spent downstream before returning to the shock 

t2^'-^. (25) 

C U2 

Consider next the other half of the cycle after the particle has crossed 
the shock from downstream to upstream until it returns to the shock. In this 
case we can define a boundary at a distance ki/ui upstream of the shock such 
that nearly all particles upstream of this boundary have never encountered the 
shock, and nearly all the particles between this boundary and the shock have 
diffused there from the shock. Then dividing the number of particles per unit 
area of shock between the shock and this boundary, ricnki/ui, by i?,cross we 
obtain the average time spent upstream before returning to the shock 

4fci 
ti « - — , 
c Ul 



(26) 



and hence the cycle time 



C \Ui U2 

The acceleration time at energy E, defined by E/{dE/dt) is then given by 



Wcle--f^ + ^). (27) 



W«-f^ + ^). (28) 



Ul \Ul U2 

We next consider the diffusion for the cases of parallel, oblique, and per- 
pendicular shocks, and estimate the maximum acceleration rate for these cases. 
The diffusion coefficients required fci and k2 are the coefficients for diffusion 
parallel to the shock normal. The diffusion coefficient along the magnetic field 
direction is some factor r] times the minimum diffusion coefficient, known as 
the Bohm diffusion coefficient, 

h = rilrgc (29) 
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where rg is the gyroradius, and 77 > 1. 

Parallel shocks are defined such that the shock normal is parallel to the 
magnetic field (_B||?ri). In this case, making the approximation that ki = k2 = 
A;|| and Bi = B2 one obtains 

II 20 f]E 

For a shock speed of Ui = 0.1c and t] — 10 one obtains an acceleration rate (in 
SI units) of 

^ « 1.5 X 10-WB. (31) 

dt acc 

For the oblique case, the angle between the magnetic field direction and 
the shock normal is different in the upstream and downstream regions, and the 
direction of the plasma fiow also changes at the shock. The diffusion coefficient 
in the direction at angle 9 to the magnetic field direction is given by 

k = kii cos^ + kj_ sin^ 9 (32) 

where k± is the diffusion coefficient perpendicular to the magnetic field. Jokipii 
13 shows that 

k^ « (33) 

provided that ry is not too large (values in the range up to 10 appear appro- 
priate). 

In the case of acceleration at perpendicular shocks, Jokipii pO| has shown 
that acceleration can be much faster than for the parallel case. For k^x — k± 
and B2 ~ 4:Bi and one obtains 



3 rjeBiuf 



For a shock speed of ui ~ 0.1c and 77 = 10 one obtains an acceleration rate (in 
SI units) of 

dE ^ 

— w OMec^B. (35) 

acc 

This discussion of shock acceleration has been of necessity brief, and has 
omitted a number of subtleties such as the finite thickness of the shock front, 
and the reader is referred to the excellent reviews cited earlier for such details. 
Nevertheless, the basic concepts have been described in sufficient detail that 
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we can consider acceleration and interactions of the highest energy cosmic rays, 
and to what energies they can be accelerated. Supernova shocks remain strong 
enough to continue accelerating cosmic rays for about 1000 years. The rate 
at which cosmic rays are accelerated is inversely proportional to the diffusion 
coefficient (faster diffusion means less time near the shock). For the maximum 
feasible acceleration rate, a typical interstellar magnetic field, and 1000 years 
for acceleration, energies of 10^^ x Z eV are possible {Z is atomic number) at 
parallel shocks and 10^^ x Z eV at perpendicular shocks. 

3 Interactions of High Energy Cosmic Rays 

Interactions of cosmic rays with radiation and magnetic fields are important 

both during acceleration when the resulting energy losses compete with energy 
gains by, for example, shock acceleration, and during propagation from the ac- 
celeration region to the observer. For ultra-high energy (UHE) cosmic rays, the 
most important processes are pion photoproduction and Bethc-Heitlcr pair pro- 
duction both on the microwave background, and synchrotron radiation. In the 
case of nuclei, photodisintegration on the microwave background is important. 
In this section I shall concentrate on photoproduction and pair production. 
The mean interaction length, Xp^ , of a proton of energy E is given by. 



where n(e) is the differential photon number density of photons of energy e, 
and a{s) is the appropriate total cross section for the process in question for a 
centre of momentum (CM) frame energy squared, s, given by 



where 9 is the angle between the directions of the proton and photon, and /3c 
is the proton's velocity. 

For pion photoproduction 




(36) 



s = mlc^ + 2eE{l - /3 cos 9) 



(37) 



/ 2 I 



;2)2 « 1.16 GeV^ 



(38) 



(39) 



2E{l + /3) AE 
For photon-proton pair-production the threshold is somewhat lower. 



s^i^ = [mpc^ + 2mec^f « 0.882 GeV^ 



(40) 
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and 

emin ~ meC^imeC^ + mpc'^)/E. (41) 

For both processes, 

smax(£, E) = my + 2sE{l + /3) « my + AeE, (42) 

and Smax(£, E) corresponds to a head-on coUision of a proton of energy E and 
a photon of energy e. 

Examination of the integrand in Equation ^ shows that the energy of the 
soft photon interacting with a proton of energy E is distributed as 



Pi^)-^Jp^n^^Ue.E)) (43) 
in the range £min < £ < oo where 

$(smax) = / a{s){s ~ my)ds. (44) 



Similarly, examination of the integrand in Equation 36 shows that the square 
of the total CM frame energy is distributed as 

<j(s){s — m'ic^) 
Pis) = f^, (45) 

in the range Smin < s < Smax- 

The Monte Carlo rejection technique can be used to sample e and s re- 
spectively from the two distributions, and Equation ^ is used to find 9. One 
then Lorentz transforms the interacting particles to the frame in which the 
interaction is treated (usually the proton rest frame), and samples momenta 
of particles produced in the interaction from the appropriate differential cross 
section by the rejection method. The energies of produced particles are then 
Lorentz transformed to the laboratory frame, and the final energy of the pro- 
ton is obtained by requiring energy conservation. In this procedure, it is not 
always possible to achieve exact conservation of both momentum and energy 
while sampling particles from inclusive differential cross sections (e.g. multiple 
pion production well above threshold) , and the momentum of the last particle 
sampled is therefore adjusted to minimize the error. 

The mean interaction lengths for both processes, Xpj{E), are obtained from 
Equation |3^ for interactions in the microwave background and are plotted as 
dashed lines in Fig. ^. Dividing by the inelasticity, k{E), one obtains the 
energy-loss distances for the two processes, 

E Xp^{E) 



dE/dx k{E) 



(46) 



11 



° 



-3 



6 



9 



12 



15 



log(£'/GeV) 



Figure 4: Moan interaction length (dashed lines) and energy- loss distance (solid lines), 
E/{dE/dx), for proton-photon pair-production and pion- production in the microwave back- 
ground (lower and higher energy curves respectively). (From Protheroe and Johnson |pl|l). 



4 Maximum Energies 

Protons and nuclei can be accelerated to much higher energies than electrons 
for a given magnetic environment. For stochastic particle acceleration by elec- 
tric fields induced by motion of magnetic fields B, the rate of energy gain by 
relativistic particles of charge Ze can be written (in SI units) 



where ^ < 1 and depends on the acceleration mechanism; a value of ^ = 0.04 
might be achieved by first order Fermi acceleration at a perpendicular shock 
with shock speed ~ 0.1c. 

^.1 Limits from Synchrotron Losses 

The rate of energy loss by synchrotron radiation of a particle of mass Amp, 
charge Ze, and energy jmc^ is 



dE 
'dt 



acc 



(47) 




(48) 
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Equating the rate of energy gain with the rate of energy loss by synchrotron 
radiation places one limit on the maximum energy achievable by electrons, 
protons and nuclei: 



El 



E' 



p 



'Z,A 



imax 



imax 



imax 



) 



-1/2 



GeV, 



GeV, 



GeV. 



(51) 



(49) 



(50) 



Other limits on the maximum energy are placed by the dimensions of the 
acceleration region and the time available for acceleration. These limits were 
obtained and discussed in some detail by Biermann & Strittmatter . 

4-2 Limits from Interactions with Radiation 

Equating the total energy loss rate for proton-photon interactions (i.e. the 
sum of pion production and Bethe-Heitler pair production) in Fig ^ to the rate 
of energy gain by acceleration gives the maximum proton energy in the absence 
of other loss processes. This is shown in Fig |^ as a function of magnetic field 
which determines the rate of energy gain through Eq. |4^. The result is shown 
by the thin curves for the maximum possible acceleration rate ^ = 1 (dashed), 
plausible acceleration at perpendicular shock ^ — 0.04 (solid), and plausible 
acceleration at parallel shock ^ = 1.5 x lO^'' (dot-dash). Also shown is the 
maximum energy determined by synchrotron losses (thick lines) for the three 
cases. As can be seen, for a perpendicular shock it is possible to accelerate 
protons to - 4 X lO^^ ggy in a - 10"^ G field. 

In the case of nuclei the situation is a little more complicated. The thresh- 
old condition for Bethe-Heitler pair production can be expressed as 



and the threshold condition for pion photoproduction can be expressed as 



Since 7 E/AmpC^, where A is the mass number, we will need to shift both 
energy-loss distance curves in Fig. ^ to higher energies by a factor of A. We 
shall also need to shift the curves up or down as discussed below. 




(52) 




(53) 
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Figure 5: Maximum proton energy as a function of magnetic field (see text). 



For Bethe-Heitler pair production the energy lost by a nucleus in each 
collision near threshold is approximately A_B « 72meC^. Hence the inelasticity 
is 

and is a factor of A lower than for protons. On the other hand, the cross section 
goes like Z^, so the overall shift is down (to lower energy- loss distance) by 
Z'^ I A. For example, for iron nuclei the energy loss distance for pair production 
is reduced by a factor 26'^/56 ~ 12.1. 

For pion production the energy lost by a nucleus in each collision near 
threshold is approximately Ai? « ^ra^^c^ ^ and so, as for pair production, the 
inelasticity is factor A lower than for protons. The cross section increases 
approximately as J^ '^ giving an overall increase in the energy loss distance for 
pion production of a factor A'^'^ k, 1.5 for iron nuclei. 

The energy loss distances for pair production and pion photoproduction, 
together with the mean free path for pion photoproduction are shown for iron 
nuclei in Fig. ^. Photodisintegration is very important and has been considered 
in detail by Tkaczyk W. et al. and Puget et al. |Q. The photodisinte- 
gration distance defined by Aj [dAj dx) taken from Puget et al. is also 
shown in Fig. ^. Since iron nuclei will be fragmented also during pion photo- 
production, the effective loss distance will be given by the photodisintegration 
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Figure 6: Mean interaction length (dashed line) and energy-loss distance (solid line), 
E/{dE/dx), for Fe-photon pion-production, and energy-loss distance for pair-production 
in the microwave background (leftmost solid line). Also shown is the photodisintcgration 
distance (thick solid curve). 

distance below ~ 3 x 10^^ GeV, and by the mean free path for pion photopro- 
duction at higher energies. The effective loss distance given in Fig. ^ is used 
together with the acceleration rate and synchrotron loss rate for iron nuclei 
to obtain the maximum energy as a function of magnetic field. This is shown 
in Fig. ^ which is analogous to Fig. ^ for protons. We see that for a perpen- 
dicular shock it is possible to accelerate iron nuclei to ~ 2 x 10"'^^ GeV in a 
^ 2 X 10^^ G field. While this is higher than for protons, iron nuclei are likely 
to get photodisintegrated into nucleons of maximum energy 400 EeV, and so 
there is not much to be gained unless the source is nearby. Of course, poten- 
tial acceleration sites need to have the appropriate combination of size (much 
larger than the gyroradius at the maximum energy), magnetic field, and shock 
velocity (or other relevant velocity), and these criteria have been discussed in 
detail by HiUas |3|. 

5 Cascading During Propagation 

There are several cascade processes which are important for UHE cosmic rays 
propagating over large distances through a radiation field: protons interact 
with photons resulting in pion production and pair production; electrons in- 
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Figure 7: Maximum iron nucleus energy as a function of magnetic field (see Fig. a). 



teract via inverse-Compton scattering and triplet pair production, and emit 
synchrotron radiation in the intergalactic magnetic field; 7-rays interact by 
pair production. Energy losses due to cosmological redshifting of high energy 
particles and 7-rays can also be important, and the cosmological redshifting of 
the background radiation fields means that energy thresholds and interaction 
lengths for the above processes also change with epoch (see e.g. Protheroe et 
al. [|6l). 

The energy density of the extragalactic background radiation is dominated 
by that from the cosmic microwave background at a temperature of 2.735 K. 
Other components of the extragalactic background radiation are discussed in 
the review of Ressel and Turner . The extragalactic radiation fields impor- 
tant for cascades initiated by UHE cosmic rays include the cosmic microwave 
background, the radio background and the infrared-optical background. The 
radio background was measured over twenty years ago [^8| , |39[| , but the fraction 
of this radio background which is truly extragalactic, and not contamination 
from our own Galaxy, is still debatable. Berezinsky was first to calculate 
the mean free path on the radio background. Recently Protheroe and Bier- 
mann jZl] have made a new calculation of the extragalactic radio background 
radiation down to kHz frequencies. The main contribution to the background 
is from normal galaxies and is uncertain due to uncertainties in their evolu- 
tion. The mean free path of photons in this radiation field as well as in the 



16 



microwave and infrared backgrounds is shown in Fig. |8| Also shown is the 
mean interaction length for muon pair-production which is negligible in com- 
parison with interactions with pair production on the radio background and 
double pair production on the microwave background. 
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Figure 8: The mean interaction length for pair production for 7-rays in the Radio Back- 
ground calculated in the present work (solid curves labelled R: upper curve - no evolution 
of normal galaxies; lower curve - pure luminosity evolution of normal galaxies) and in the 
radio background of Clark (dotted line). Also shown are the mean interaction length for 
pair production in the microwave background (2.7K), the infrared and optical background 
(IR), and miiori pair production and double pair production (4e) in the microwave 

background pi| . (From Protheroe and Biermann jll|). 

Inverse Compton interactions of high energy electrons and triplet pair 
production can be modelled by the Monte Carlo technique ji^ |4^], and the 
mean interaction lengths and energy-loss distances for these processes are given 
in Fig. p. Synchrotron losses must also be included in calculations and the 
energy-loss distance has been added to Fig. ^ for various magnetic fields. 

5.1 Practical Aspects of the Cascade 

Where possible, to take account of the exact energy dependences of cross- 
sections, one can use the Monte Carlo method. However, direct application 
of Monte Carlo techniques to cascades dominated by the physical processes 
described above over cosmological distances takes excessive computing time. 
Another approach based on the matrix multiplication method has been de- 
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Figure 9: The mean interaction length (dashed line) and energy- loss distance (solid line), 
E/{dE/dx), for electron-photon triplet pair production (TPP) and inverse-Compton scatter- 
ing (IC) in the microwave background. The energy-loss distance for synchrotron radiation is 
also shown (dotted lines) for intergalactic magnetic fields of 10~® (bottom), 10"^", 10~^^, 
and 10~^^ gauss (top). (From Prothcroe and Johnson |pi[). 



scribed by Protheroe and developed in later papers A Monte 

Carlo program is used to calculate the yields of secondary particles due to 
interactions with radiation, and spectra of produced pions are decayed (e.g. 
using routines in SIBYLL [^Tf) to give yields of gamma-rays, electrons and neu- 
trinos. The yields are then used to build up transfer matrices which describe 
the change in the spectra of particles produced after propagating through the 
radiation fields for a distance 5x. Manipulation of the transfer matrices as 
described below enables one to calculate the spectra of particles resulting from 
propagation over arbitrarily large distances. 



5.2 Matrix Method 

In the work of Protheroe and Johnson |3l| ] , fixed logarithmic energy bins were 
used, and the energy spectra of particles of type a {a — ^, e,p, n, v^, f'e, v^, i>^) 
at distance x in the cascade are represented by vectors F^{x) which give the 
total number of particles of type a in the jth energy bin at distance x. Transfer 
matrices, T'^^{5x), give the number of particles of type (3 in the bin j which 
result at a distance 5x after a particle of type a and energy in bin i initiates a 
cascade. Then, given the spectra of particles at distance x one can obtain the 
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spectra at distance {x + Sx) 

180 

F^^{x + Sx)^Y.T.^^f(^^)P"('^) (55) 

a i=j 

where F°'{x) are the input spectra (number in the ith energy bin) of species 
a. 

We could also write this as 

[F{x + Sx)] = [T{Sx)] [F{x)] (56) 

where 





pi 




^71 


1^67 


. . . " 


















[F] = 


pP 


, [T] = 








(57) 



The transfer matrices depend on particle yields, , which are defined 
as the probability of producing a particle of type f3 in the energy bin j when 
a primary particle of type a with energy in bin i undergoes an interaction. 
To calculate Y^'' a Monte Carlo simulation can be used (see Protheroe and 
Johnson for details). 

5.3 Matrix Doubling 

From Fig. || we see that the smallest effective interaction length is that for 
synchrotron losses by electrons at high energies. We require Sx be much smaller 
than this distance which is of the order of parsecs for the highest magnetic field 
considered. To follow the cascade for a distance corresponding to a redshift of 
z ~ 9, and to complete the calculation of the cascade using repeated application 
of the transfer matrices would require ~ 10^^ steps. This is clearly impractical, 
and one must use the more sophisticated approach described below. 

The matrix method and matrix doubling technique have been used for 
many years in radiative transfer problems [^^9[|. The method described by 
Protheroe and Stanev ji^ is summarized below. Once the transfer matrices 
have been calculated for a distance Sx^ the transfer matrix for a distance 2Sx 
is simply given by applying the transfer matrices twice, i.e. 

[T{2Sx)] = [T{Sx)]^. (58) 

In practice, it is necessary to use high-precision during computation (e.g. 
double-precision in FORTRAN), and to ensure that energy conservation is 
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preserved after each doubling. The new matrices may then be used to calcu- 
late the transfer matrices for distance ASx, 86x, and so on. A distance 2"fe 
only requires the application of this 'matrix doubling' n times. The spectrum 
of electrons and photons after a large distance Ax is then given by 

[F(a; + Ax)] = [T(Aa;)][F(a;)] (59) 

where [F(a;)] represents the input spectra, and Ax = 2^Sx. In this way, cas- 
cades over long distances can be modelled quickly and efficiently. 

6 The Origin of Cosmic Rays between 100 TeV and 300 EeV 

The subject of the origin of cosmic rays at these energies has been reviewed 
[^,^, and one of the very few plausible acceleration sites may be associated 
with the radio lobes of powerful radio galaxies, either in the hot spots or 
possibly the cocoon or jet One-shot processes comprise another possible 
class of sources p3| , |5^ . 

Acceleration at the termination shock of the galactic wind from our Galaxy 
has been also been suggested Q , but due to the lack of any statistically signif- 
icant anisotropy associated with the Galaxy is unlikely to be the explanation. 
However, a very recent re-evaluation of the world data set of cosmic rays has 
shown that there is a correlation of the arrival directions of cosmic rays above 
40 EeV with the supergalactic plane (56) , lending support to an extragalactic 
origin above this energy, and in particular to models where "local" sources 
(< 100 Mpc) would appear to cluster near the supergalactic plane (e.g. pow- 
erful radio galaxies as in the model of Rachen and Biermann ||5l|] . 

Rachen and Biermann have demonstrated that cosmic ray acceleration 
in Fanaroff-Riley Class II radio galaxies can fit the observed spectral shape 
and the normalization at 10 - 100 EeV to within a factor of less than 10. The 
predicted spectrum below this energy also fits the proton spectrum inferred 
from Fly's Eye data Protheroe and Johnson have repeated Rachen 
and Biermann's calculation to calculate the flux of diffuse neutrinos and gamma 
rays which would accompany the UHE cosmic rays, and their result is shown 
in Fig. 1^. 

6.1 Observability of Ultra High Energy Gamma Rays 

Above 100 EeV the interaction properties of gamma rays in the terrestrial 
environment are very uncertain. Two effects may play a significant role: in- 
teraction with the geomagnetic field, and the Landau-Pomeranchuk-Migdal 
(LPM) effect |5|,|60). 
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Figure 10: Cosmic ray proton intensity multiplied by in the model of Raclien and Bier- 
mann for Hq = 75 km s~^ Mpc~^ with proton injection up to 3 X 10^^ GeV (solid line). 
Also shown are intensities of neutrinos (dotted lineSj.i'^, P^, fe, i'e from top to bottom), and 
photons (long dashed lines). Data are from Stanev [p8|; large crosses at EeV energies are an 
estimate of the proton contribution to the total intensity based on Fly's Eye observations. 
(From Protheroe and Johnson pi|). 



Energetic gamma rays entering the atmosphere will be subject to the LPM 
effect (the suppression of electromagnetic cross-sections at high energy) which 
becomes very important at ultra-high energies. The radiation length changes 
as {E/Ei^puy^^, where E'lpm — 6.15 x lO^fcm GeV, and 4m is the standard 
Bethe-Heitler radiation length in cm Protheroe and Stanev |^ found 

that average shower maximum will be reached below sea level for energies 
5 X 10" GeV, 8 x 10" GeV, and 1.3 x lO^^ Q^y for gamma-rays entering 
the atmosphere at cos 6* = 1, 0.75, and 0.5 respectively. Such showers would 
be very difficult to reconstruct by experiments such as Fly's Eye and at best 
would be assigned a lower energy. 

Before entering the Earth's atmosphere gamma rays and electrons are 
likely to interact on the geomagnetic field (see Erber for a review of the 
theoretical and experimental understanding of the interactions). In such a case 
the gamma rays propagating perpendicular to the geomagnetic field lines would 
cascade in the geomagnetic field, i.e. pair production followed by synchrotron 
radiation. The cascade process would degrade the gamma ray energies to some 
extent (depending on pitch angle), and the atmospheric cascade would then 
be generated by a bunch of gamma-rays of lower energy. Aharonian et al. []64| 
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have considered this possibihty and conclude that this bunch would appear 
as one air shower made up of the superposition of many air showers of lower 
energy, where the LPM effect is negligible; the air shower having the energy of 
the initial gamma-ray outside the geomagnetic field. If this is the case, then 
gamma rays above 300 EeV would be observable by Fly's Eye, etc. There is 
however, some uncertainty as to whether pair production will take place in the 
geomagnetic field. This depends on whether the geomagnetic field spatial di- 
mension is larger than the formation length of the electron pair, i.e. the length 
required to achieve a separation between the two electrons that is greater than 
the classical radius of the electron (see also Stanev and Vankov ||65|). 

7 Exotic Origin Models 

We now discuss the possibility that the highest energy cosmic rays are not 
single nucleons. Obvious candidates are heavier nuclei (e.g. Fe), 7-rays, and 
neutrinos. In general it is even more difficult to propagate nuclei than protons, 
because of the additional photonuclear disintegration which occurs ||3^,|3^,|6^. 
The possibility that the 300 EeV event is a 7-ray has been discussed recently 
and, although not completely ruled out, the air shower development profile 
seems inconsistent with a 7-ray primary. Weakly interacting particles such as 
neutrinos will have no difficulty in propagating over extragalactic distances, of 
course. This possibility has been considered, and generally discounted p^,|67|, 
mainly because of the relative unlikelihood of a neutrino interacting in the 
atmosphere, and the necessarily great increase in the luminosity required of 
cosmic sources. Magnetic monopoles accelerated by magnetic fields in our 
Galaxy have also been suggested and can not be ruled out as the highest 
energy events until the expected air shower development of a monopole-induced 
shower is worked out. 

7.1 Topological Defect Origin 

One exciting possible explanation of the highest energy cosmic rays is the 
topological defect (TD) scenario ||69|-|7l|], where the observed cosmic rays are a 
result of top-down cascading, from somewhat below GUT scale energy, ~ 10^^ 
GeV, to 10^^ GeV and lower energies. Generally, these models put out much of 
the energy in a very flat spectrum of photons and electrons extending up to the 
mass of "X-particles" emitted which may be lower than 10^^ GeV, depending 
on the theory. Approximating this spectrum by a monoenergetic injection of 
photons of energy 10^^ GeV, Protheroe & Johnson |^l|] showed that spectra 
from single TD sources can not explain the 3 x 10^^ GeV events. 
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The main problem with topological defect models is the wide range of 
model parameters in which this scenario could, in principle, be applied. This 
is very different from astrophysical scenarios [pT| where nearly all of the prop- 
erties of the astrophysical objects are restricted from observations (luminosities 
at different wavelengths, magnetic field, etc.). Parameters of TD scenarios in- 
clude: mass of the X-particle (maximum injection energy for X-particle decay 
products), energy spectra and final state composition of the decay products, 
and cosmological evolution of the topological defect injection rate |72,73|. 

To compare the TD decay products with the cosmic ray data one has to 
study the change of the spectra during propagation |^,|7^. This problem is 
more severe than for the case of "bottom-up" scenarios because most of the 
energy from X-particle decay, and subsequent decays, emerges in electrons, 
photons and neutrinos, with only about 3% in nucleons. The electrons and 
photons initiate electromagnetic cascades in the extragalactic radiation fields 
and magnetic field, resulting in a complicated spectrum of electrons and pho- 
tons which is very sensitive to the radiation and magnetic environment. For 
example, recently the HEGRA group have placed an upper limit on the 
ratio of gamma-rays to cosmic rays of ^ 10~^ at 10^ GeV, and based on a TD 
model calculation which neglected the IR background and gave a higher 
ratio argued that TD models were ruled out. However, inclusion of the IR 
reduces the 10^ GeV gamma-ray intensity to below the HEGRA limit. 

Unfortunately, many of the relevant parameters of extragalactic space (av- 
erage magnetic field, strength of the radio and IR/optical backgrounds) are not 
well known. Protheroe and Johnson have considered one set of parame- 
ters {mxc^ — 10^^ GeV, constant injection per co-moving volume, B — 10~^ 
gauss) and rule out TD as the origin of the 3 x 10^^ GeV events. More recently, 
Lee Q and Sigl, Lee and Coppi |Q, adopting a probably unphysical lower 
X-particle mass mx(? — 10^^ GeV (the X-particle mass is expected to be near 
the unification mass which is 10i6-0±o.3 ^ iQ^gj. 

magnetic field 

claim the TD scenario is not ruled out. However, Protheroe and Stanev 
argue that these X-particles masses are ruled out as well. The details are not 
simple, and include considering whether or not UHE gamma rays in the cas- 
cade are observable or not by air shower arrays, but in either case TD models 
appear to be ruled out as the origin of the 300 Ee V events. 



8 Conclusion 

While the origin of the highest energy cosmic rays remains uncertain, there 
appears to be no necessity to invoke exotic models. Shock acceleration, which 
is believed to be responsible for the cosmic rays up to at least 100 GeV, is a 
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well-understood mechanism and there is evidence of shock acceleration taking 
place in radio galaxies , and the conditions there are favourable for acceler- 
ation to at least 300 EeV. Such a model, in which cosmic rays are accelerated 
in Fanaroff-Riley Class II radio galaxies, can readily account for the flat com- 
ponent of cosmic rays which dominates the spectrum above ~ 10 EeV p7| . 
Indeed, one of the brightest FR II galaxies, 3C 134, is a candidate source for 
the 300 EeV Fly's Eye event (P.L. Biermann, personal communication). 

Whatever the source of the highest energy cosmic rays, because of their 
interactions with the radiation and magnetic fields in the universe, the cosmic 
rays reaching Earth will have spectra, composition and arrival directions af- 
fected by propagation. Many of the important parameters, e.g. extragalactic 
magnetic fields and radiation fields, are uncertain, and more work is needed to 
obtain better estimates of these and thereby help unravel the puzzle. When 
better information is available, and good statistics on the arrival directions, 
energy spectra and composition, as well as the intensity of the diffuse back- 
ground of very high energy gamma rays (partly produced in cascades initiated 
by cosmic ray interactions), we will be better placed to understand the origin 
of the highest energy particles occurring in nature. The Auger Project, an 
international collaboration to build two UHE cosmic ray detectors, one in the 
United States and one in Argentina, each having a collecting area in excess of 
1000 km^ (see the lecture by J.M. Matthews in this volume), will go a long 
way to help answer these questions. 
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